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Abstract. The intertwining technique has been widely used to study the Schrodinger 
equation and to generate new Hamiltonians with known spectra. This technique can be 
adapted to find the bound states of certain Dirac Hamiltonians. In this paper the system 
to be solved is a relativistic particle placed in a magnetic field with cylindrical symmetry 
whose intensity decreases as the distance to the symmetry axis grows and its field lines are 
parallel to the x — y plane. It will be shown that the Hamiltonian under study turns out to 
be shape invariant. 
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1 Introduction 

The intertwining technique, also called Supersymmetric Quantum Mechanics (SUSY QM), is 
a widespread method used to generate exactly solvable Hamiltonians departing from a given 
initial one and can be employed as well to solve a certain set of Hamiltonians in a closed way, 
among other applications. In the simplest case (1-SUSY QM) the new potentials have similar 
spectra as the original one, namely, they might differ at most in the ground state energy. 
Examples of potentials generated by this technique are those which arise when adding a bound 
state to the free particle Hamiltonian (hyperbolic Poschl-Teller) [14] or the Abraham-Moses- 
Mielnik potentials which are isospectral to the harmonic oscillator [1, 13, 15]. This method has 
been also applied successfully to the radial part of the hydrogen atom potential [1, 7, 13, 18], 
the trigonometric Poschl-Teller potentials [3], among many others. 

To apply the technique [8] we start from two one-dimensional Schrodinger Hamiltonians 



where Hq is known. Now let us suppose the existence of a differential operator A[ which satisfies 



Since the operator A[ is of first order, the technique is known as 1-SUSY QM and the func- 
tion W\(x) as the superpotential. It is also said that the potentials Vq{x) and V\(x), whose 
Hamiltonians are intertwined by the operator A\, are supersymmetric partners. 



1 d 2 



i = 0, 1 





(1) 
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In order to satisfy equation (1) V\{x) and W\(x) must obey 

V 1 {x) = V (x)-W[(x), W[(x) + W?(x) = 2(V -ei), (2) 

where e\ is a real integration constant called factorization energy. From the previous equations 
we can see that if Vq{x) is given and Wi(x) is found, the supersymmetric partner V\{x) is 
completely determined. Furthermore, equation (1) ensures that if ip n is an eigenfunction of Hq 
with eigenvalue E n then A\ip n will be an eigenfunction of H\ with the same eigenvalue. Note 
that the operators A\ and {A\)^ = A\ factorize the Hamiltonian as follows 

H = AiA[ + ex, H x = A\A x + e u (3) 

where 

By taking the squared norm of the vectors A\ip n we have ||A[^ n || 2 = {A\ip n , A\ip n ) = 
(ip n , AiA\ifj n ) = E n — e\ > Vn which implies that e\ < Eq, where Eq is the ground state 
energy of Hq. One could ask now if {A\iP n , n = 0,1,2,...} is a complete orthogonal set. In 
order to answer this question, assume the existence of a vector ip ei orthogonal to each vector of 
the previous set, then 

{lp ei , A\lP n ) = (AnP tl ,?/; n ) = Vn Ailpe, = 0, 

since {ip n , n = 0, 1, 2, . . . } is a complete orthogonal set. The first-order differential equation 
A\ip ei = can be solved immediately 

/•X 

tp ei oc exp - / Wi(y)dy 
. J o 

Notice that ip ei satisfies 
Hiil> €1 = eitl> ei . 

Thus, depending on the square integrability of this vector, and the value of ei, three possibilities 
arise 

• The function ip ei with ei < Eq belongs to the Hilbert space H. Thus, {ip ei , A\ip n , n = 
0, 1, 2, . . . } is a complete orthogonal set, and from equations (1), (3) the spectrum of H\ 
is given by Sp[#i] = {ei, E n , n = 0, 1, 2, . . . }. 

• ip ei ^ 71 with ei < E\. In this case {A\ip n , n = 0, 1, 2, . . . } is a complete orthogonal set 
and thus Sp[i/i] = Sp[i/o]- 

• When tp ei ^ H, t\ = Eq, the set {A[ip n , n = 1, 2, 3, . . . } is complete and thus Sp[i^i] = 
{E n , n = 1,2,3,...}. 

Restricting ourselves to this last case, it can be verified that Wi(x) = iP'q/iPq fulfills equa- 
tion (2) and applying successively this technique we can generate a hierarchy of Hamiltonians, 
where Spfi^o] D Sp[f/i] D Spf^] ^ • • • • This sequence is either finite or infinite if the number 
of bound states of Hq is finite or infinite respectively [21]. 

Up to this point we have assumed that starting from a solvable Hamiltonian Hq we can ge- 
nerate a hierarchy of Hamiltonians {Hi, i = 0,1, . . .}. However, from the equation adjoint to 
equation (1) we can see that beginning from an eigenvector of H{ we can construct an eigenvector 
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of ifj-i through the action of the operator A{. Thus, if we had known enough eigenvectors of 
the Hamiltonians of the hierarchy, for example all the ground states, it would be possible to 
build all bound states of Hq by applying the operators Ai over them. 

It is convenient to recall now the concept of shape invariance. If two SUSY partner potentials 
Vi,2(x;ai) satisfy the condition 

V 2 (x; 01) = Vi(x; a 2 ) + R(ai), 

where a\ is a set of parameters, a 2 is a function of a\ and the remainder R{ai) is independent 
of x, then V\(x\ a±) and V 2 (x; a{) are said to be shape invariant [4]. 

If the potentials of a hierarchy of Hamiltonians are shape invariant and the ground state of 
one of them is found, in principle, all the ground states can be derived. In this way we can find 
the eigenfunctions of the first Hamiltonian. The harmonic oscillator and the radial effective 
potential of the hydrogen atom are examples of shape invariant potentials that can be solved 
through this procedure. 

It is noteworthy that there are papers in which through the SUSY technique the Dirac 
equation for different systems has been solved [2, 5, 6, 12, 16, 17, 20] or analyzed [11]. However, 
the differences with respect to the approach we will use here will be significant. 

In Section 2 we will employ the 1-SUSY QM in order to solve the stationary Schrodinger 
equation for a charged particle placed in a magnetic field generated by the vector potential 
A(p, <ft, z) = ^e z , where k is a constant characterizing the field strength, c is the speed of light, 
e is the charge of the particle, and p is the radial variable in cylindrical coordinates. The resulting 
field has cylindrical symmetry, its intensity decreases with the distance to the symmetry axis 
and its field lines are parallel to the x — y plane. Making use of the basic ideas to solve the shape 
invariant potentials through the 1-SUSY QM technique, we will work out the same problem in 
Section 3 in the relativistic regime by solving the associated Dirac equation. In the last section 
we will present our conclusions. 



2 Nonrelativistic quantum approach 

The classical Hamiltonian of a particle with charge e and mass m in a magnetic field generated 
by the vector potential A(p, cf>, z) = ^e 2 is given by 

2m 2m 2m \ p J 

The corresponding magnetic field B = V x A could be produced in a coaxial transmission line, 
with the inner and outer conductors carrying currents I a and —I a T~/(r + r) respectively, where 
r Q is the minor radius of the outer conductor and r is its thickness. If the current density in 
the second conductor is ||J|| = r Q I a / '(27T/5 3 ) then such a magnetic field will be generated in the 
material [9, 19]. 

In order to address the quantum treatment the classical observables have to be promoted to 
the corresponding quantum operators. The quantum Hamiltonian is then 

2m mp \ oz J 2mp z 

where V 2 is the Laplacian operator. 

It can be seen that the Hamiltonian commutes with the operators of partial derivative with 
respect to z and 0, 
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-b 2 / [2(a + 4) 2 ] 
-b 2 / [2(a + 3) 2 ] 
-b 2 / [2(a + 2) 2 ] 
-b 2 / [2(a + l) 2 ] 



H Hi H 2 H3 

Ai A 2 A 3 



G M = A x A 2 A 3 <t> H A 2 A 3 (f> ei 



G-21 = A^fc, A 2 (j> t 



Got = </>ei 



A 



<!>,: 



Figure 1. A hierarchy of Hamiltonians built up departing from Hq. If we know the ground state of each 
Hamiltonian and the intertwining operators, we can know the bound states of all Hamiltonians. Note 
that the eigenvalues are not indeed equidistant. 



This suggests us the following ansatz for the solutions of the stationary Schrodinger equation 

i;(p,4>,z) = e^ z+ ^/ h p- 1 / 2 G(p), (5) 

where p z and £ are respectively the eigenvalues of the momentum operator along z, P z = 
—ihd/dz, and the z component of the angular momentum, L z = —ihd/d<p. 

Through this ansatz we can separate variables for the stationary Schrodinger equation, Hip = 
Eip, leading us to the differential equation for G(p) 



1 d 2 (A/ft) 2 - 1/4 Pz k 

2 dp 2 2p 2 h 2 p 



GM-S(*-£|GM, 



where A 2 = i 2 + k 2 . To simplify notation we can express the previous equation as 



1 d 2 a(a + 1) 



2dp 2 



+ 



2p 2 



G(p) = dG(p), 



(6) 



with 



a(a + 1) 



h 2 



h 2 ' 



d 



mE p 2 



In this work we restrict ourselves to the case p z k > 0, which is the one with bound states. 
Equation (6) can be identified as the radial equation of the hydrogen atom. To solve this equation 
we propose the existence of a family of operators An that intertwine the Hamiltonians H n 
and H n+ \ in the way 



H n +iA n+1 — A n+1 H n , 



where 



H„ 



Id 2 
2 dp 



, 2 + V n (p) 



1 d 2 (a + n)(a + n + 1) 
2df? + 2~i 2 



b 
P 



and 
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From equations (2) we have 
xxr ( \ a + n b 

W n (p) : 



p a + n 1 " 2{a + n) 2 

Looking for the function annihilated by A n , which due to equation (3) is an eigenfunction 
of -ff n _i with eigenvalue e n , 

it turns out that the ground state of H n is given by 

As expected, one can verify that 

W n+1 (p) = & n+1 (p)/cp tn+1 (p). (7) 

It is enough to know the ground state and its eigenvalue for any Hamiltonian of the hierarchy 
in order to find the complete solution of Hq (see Fig. 1). The spectrum is given by 

Sp[H } = \- / n = 0,l,2,.. 

[ 2(a + n + \) £ 

and its eigenfunctions by 

Gw = (j)e 1 , Gu = Ai<f> e2 , G 2 e = AiA 2 <f)e3, G 3i = AiA 2 A 3 (f> e4 , 

where the index n indicates the energy level of Hq and the index I reminds us that the radial 
Hamiltonian depends on the angular momentum. In Fig. 2 the first three eigenfunctions of Hq 
can be seen (black continuous, dashed and dotted lines) placed at its corresponding energy level 
and the potential Vq(p) is as well drawn (gray line). 

Returning to the original problem, i.e. the eigenvalue equation for the operator of equation (4), 
we have that the spectrum is given by 



Sp [H] 



2m 



1 



, rc = 0,l,2,... 



h 2 (X/h + n + l/2) 2 _ 
and its eigenfunctions by 

where C n e Pz is a normalization constant (there is not sum convention). 

3 Relativist ic quantum approach 

The stationary Dirac equation of a free particle with mass m and spin 1/2 is 

n D $> = [ca-P + f3mc 2 ]^ = Ety, (8) 

where P is the momentum operator, ^ is a four-component spinor and m and (3 are 4x4 
matrices given by 

Wi J ' H I -co 
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Vo( P ) 

-b 2 / [2(a + 3) 2 ] 
-b 2 / [2(a + 2) 2 ] 

-b 2 / [2(a + l) 2 ] 




Figure 2. The potential Vo(p) (gray curve) and its first three eigenfunctions, Go£ (black continuous 
line), Gu (dashed line) and Gn (dotted line), with a — 1.5 and b = 0.5 in units of 1/p. 



being <jq the 2x2 identity matrix and <7j the Pauli matrices. The 4x4 matrix operators are 
written in boldface in order to be distinguished from the 2x2 matrix operators. In our case 
the interaction with the magnetic field derived from the vector potential A = ^e z is described 

by the minimal coupling rule P — > P — ^A. In cylindrical coordinates the resulting stationary 
Dirac equation is 

H fl $ = <^ -i/icD(0)ai- D(</>)a 2 — - ihca 3 - a 3 + /3mc 2 } V = EV, (9) 

[ op p oq> oz p J 

where D(0) = Diag [e _i< ^, e 1 ^, e~ l< ^ , e l< ^~\ is a diagonal matrix. This interacting Hamiltonian com- 
mutes with the momentum operator, and with the total angular momentum in the z-direction, 

p 2 = -ihd z i, 3 Z = -ihd+i + Js 3 , = r* 

Then we will look for a solution to equation (9) that is also an eigenfunction of these two 
operators with corresponding eigenvalues p z and I respectively (see equation (5)), having the 
form 

*M,z) = e*'VV'(*- E »/ 2 WV 1/a Gi>(p)- 
The equation that the radial function Gd(p) must fulfill is 



d £ ( k p z \ mc 

-iai— + — a 2 -It r- «3 + — p 

dp hp \hp h J a 



G D (p) = §Gd{ P ). 



The operator between brackets is an effective Hamiltonian that will be called Hp. It is useful to 
perform an unitary transformation in order to leave all the dependence of p on a single matrix, 

H = UjHpUi, Ui = e"^ El/2 = cos(0/2)l - i sm(0/2)Ei, (10) 

with tan (9 = —k/£. In the rotated frame the Hamiltonian is 

d / A p z k\ p z £ mc 

H ° = ~ iai dp + [np--nx) a2+ nx as + ~h 

where in this context once again A 2 = £ 2 + k 2 . This Hamiltonian has a special structure that 
can be better appreciated if we write it as follows 

/ (mc/h)a h \ , d f X p z k\ p z £ 

V ho -{mc hjaoj dp \hp h\ J h\ 
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being ho a 2 x 2 matrix operator. In order to simplify notation we will write 



. d 

dp 



with 



a = X/H, 



a _ bo 
p a 

Pzk/fr 2 



<?2 +d cr 3 = -20i — + v (a,bo,d ;p), 
dp 



do = Pz^/h\. 



To solve the eigenvalue equation for Ho with a method similar to that used in the nonrelativistic 
approach, we propose the intertwining relationship 



TJ a t —At I I 

n "+l A n+l — ^n+l"-r 



(11) 



with the intertwining operators A^ +1 and the sequence of Hamiltonians H n +i having the form 



A f 

A n+1 



B 



n+l 





B 



H, 



n+l, 



(mc/h)ao 



-{mc/K)ao J ' 



(12) 



where h n (a n ,b n ,d n ; p) = ho(a + n,b n ,d n ; p) with parameters b n and d n to be determined 
and Br n , 1 is a 2 x 2 operator that intertwines /i n +i with /i n . The structure of A^ +1 above 
proposed is the simplest choice. A more general form of the intertwining operators will be 
presented elsewhere. Similar 2x2 intertwining operators were considered in [10]. 

In the same way as in the nonrelativistic quantum case (see equation (1)), the first order 
intertwining operator has the following structure 



A f 

A n+1 



where W^, ^(p) is a variable matrix. In the same way sLj reads 



B 



t 

n+l 



dp 



with F n+ \(p) a 2 x 2 matrix. 

Solving the intertwining relation of equation (11) we find 



, , 2 _ ,2 n ( 2a + n )b 

o, a n — «q + n2 / n i „\2 



5 



n+l 



d 2(a + n) + l 
+ 2 



a 2 (a + n) 2 
1 



1 



+ 



a + n)(a + n + 1) 
03 • 



00 



%+i 



C?n 



(icri - cr 2 ) 



p (a + n)(a + re + 1)_ 

Then the intertwining operators Ajj +1 are directly obtained by substitution in equation (12). 

Now, we look for the vector wavefunctions annihilated by A^ +1 . Taking advantage of the 
block diagonal structure of this operator, first we find the two-component functions annihilated 
by These are given by the following two independent vector functions 



Xr, 



P 



a+n —bp/(a+n) 



. (a + n) 2 (a + re + l) 2 (d n+1 - d n ) 



(a + n){a + n + 1) 



pi+n e - bp I (a+n+ 1 ) 
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We can built four-component functions annihilated by , 1 such that at the same time they are 
also eigenvectors of H n . We find four types of such vectors and the corresponding eigenvalues: 



Xr 



y \J (mc/h) 2 + d n + mc/h 
( 



Xn 



Xr 



\ 



'Xr, 



\ 



y yj (mc/h) 2 + d 2 — mc/h j 

dn+1 f. 

i = sn 

\ yj (mc/h) 2 + d 2 n+l + mc/h J 

in \ 



T-n+1 



; sn 

\yj (mc/h) 2 +d 2 n+l - mc/h J 



yj(mc/h) 2 + d 2 ; 



V(mc/h) 2 +d 2 ; 



^(mc/h) 2 + dl +1 ; 



^(mc/h) 2 +d 2 n+1 . 



We will briefly comment here on some properties of these results, a more complete discussion 
will be given elsewhere. 

• Superpotential. Consider now the matrix which is constructed by placing the 
vectors 4> un in its columns, it can be verified, in analogy with equation (7), that 

W 1+1(P) = -K+l(p) E nll(p)- 

• Spectrum. From the intertwining relationship, equation (11), it is shown that the spectrum 
of Hq is given by 



Sp [H c 



mc\ 2 

~h) 



2 n(2a + n)b 2 
) + + a 2 (a + n) 2 ' 



71 = 0,1,.... 



Eigenf unctions. The eigenfunctions of the initial Hamiltonian are computed in the usual 
way. We have four types of eigenfunctions: 



&a0£ = <AaO, 


®all 


= A 1( 


t>al, 


$a2f 


= AiA 2 ( 


i>a2 


$b0£ = <j>bO, 


$bU 


= Ai<) 


^61, 


$b2£ 


= AiA 2 <; 


k2, 


^c0£ = 4>c0, 


*cW 


= Ai<; 




$ C 2£ 


= AiA 2 <; 


i>c2 




*e£U 


= Ak 




$d2£ 


= AiA 2 ( 


t>d2 



where we add the subindex I to remind that the Hamiltonian Ho depends on I. 
Then the eigenvectors of the Dirac equation, equation (8), for our system are 
VvrUpMh*) = C^ P2 e^^ ,! 'e i ^-5 s 3)0/fi Ul/O -i/2 ( j ) ^ (/o)) 

where Ui is given by equation (10), v = a,b,c,d, n = 0, 1, 2, . . . and C un g Pz are normalization 
constants (there is not sum convection). The spectrum of H/j is ch times the one of Hq, thus 



Sp [H 



D 



±mc 2 i 1+ '" 



P 



Ik 2 



m 2 c 2 h 2 m 2 c 2 (\/h + n) 2 



71 = 0,1,2,.... 
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P 



(a) (b) 

Figure 3. Probability densities for six eigenvectors of Hq: in (a) the first three of the family & an e, and 
in (b) the first three of the family $ c „£, the parameters used were a = 1; b = 2 and do = 1 with units 
of 1/p; m = 0.1 and c = H = 1. 

In Fig. 3 we can see the probability densities of six eigenvectors of Ho- In (a) we have the first 
three with subindex a, and in (b) the corresponding but with subindex c. Note that & a u an d & c oi 
have the same eigenvalue but the behavior of the probability density is quite different, the same 
happens with & a 2£ and $ c i£, and so on. This degeneracy, which does not appear in the nonrela- 
tivistic approach, is due to the spin degree of freedom and will be analyzed in detail elsewhere. 

4 Conclusions 

In this work we have adapted the intertwining technique to solve exactly the Dirac equation asso- 
ciated to a charged particle of spin 1/2 immersed in a magnetic field with cylindrical symmetry 
generated by the vector potential A = ^e z . We first addressed the problem in the nonrelativistic 
regime, i.e., the Schrodinger equation through the standard intertwining technique. Afterwards 
we set up the corresponding Dirac equation and we proposed, as in the nonrelativistic approach, 
a hierarchy of shape invariant Hamiltonians intertwined by some operators to be determined. 
These operators afterwards were found and using them the ground states of each Hamiltonian 
were built. Applying these operators onto the ground states all the bound states were obtained 
as well as their respective eigenvalues of the original Dirac equation. As far as we know these 
solutions have not been reported before. The analogies between the method for solving the Dirac 
equation and the standard intertwining technique for the Schrodinger equation were recurrently 
employed throughout the entire procedure. 
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